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Nuclear effects in neutrino-nucleus reactions simulated by means of the NUX+FLUKA Monte Carlo generator 
are compared with the theoretical predictions of the Marteau model. Pion absorption in NUX+FLUKA and 
non-pionic A decays in the Marteau model differ by about 30%. The fraction of pions produced due to the 
re-interactions after primary quasi-elastic vertex is in the NUX+FLUKA scheme much higher then provided by 
the Marteau model. 



1. INTRODUCTION 

The aim of our investigation was to compare 
predictions for neutrino interactions provided 
by NUX+FLUKA Monte Carlo scheme [T] with 
those of the Marteau model PJ. We were inter- 
ested mostly in nuclear effects in pion production 
but for completeness (and for the sake of checking 
normalization factors) we have also presented the 
results for the quasi-elastic reaction. 

The motivation for this study can be summa- 
rized in the following points: 

i) In recent years there has been growing in- 
terest in the studies of neutrino interactions at 
energies of a few GeV where the dominant con- 
tributions to the cross section come from quasi- 
elastic and single pion production channels 

ii) It is well known that the NUX+FLUKA MC 
scheme designed originally to describe neutrino 
interactions at higher energies is not satisfactory 
in a few GeV neutrino energy region 4 . It does 
not include separate resonance contribution im- 
plemented in alternative MC codes after Rein- 
Sehgal model jSj. 

iii) Developments in quark-hadron duality sug- 
gest that perhaps it is not necessary to include 
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many resonances apart from the A excitation 

Ed 

iv) The Marteau model describing together 
quasi-elastic and A production reactions on nu- 
clear targets is an interesting candidate for im- 
plementation in MC codes. Its original version 
includes a contribution from 2 particles - 2 holes 
excitations. The model was a subject of further 
investigation [5] and it is important to understand 
if it is useful for practical applications. 

v) On the purely theoretical side there are in- 
teresting questions concerning the way in which 
nuclear effects beyond Fermi gas model are taken 
into account. One way is to describe them nu- 
merically e.g. by FLUKA. An alternative is to 
perform sophisticated theoretical computations 
like RPA. It is important to compare predictions 
of these two approaches. 

In this paper we confine ourselves to the CC 
reactions of muon neutrinos. The NUX+FLUKA 
Monte Carlo predictions were made by using the 
PRET driver. For the total cross sections plots 
a sample of 10 5 events was produced for each 
neutrino energy value. The events were classified 
based on some primary vertex characteristics as 
well as the number of pions in the final state. The 
overall normalization factor for the cross section 
was taken from the PRET's output file prel.out. 
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Figure 1. The total cross section for the CC quasi- 
elastic Vy, scattering on free nucleon. The experi- 
mental points are taken from 0. 



Figure 2. The total cross section for the CC quasi- 
elastic Vy scattering on 16 O without Final State 
Interactions. 



The differential cross section, i.e. the hadronic 
invariant mass distribution in the vp — > /i _ p7r + 
channel, was obtained from a sample of 10 6 in- 
clastic events with the neutrino energy fixed at 
lGeV. 

All the numerical predictions of the Marteau 
model are based on [S]- 

2. QUASI-ELASTIC REACTIONS 

The total cross sections for CC quasi-elastic 
scattering on free nucleons are compared in Fig. 1. 
We find good agreement between the two plots. 
The small difference may be attributed to slightly 
different parameter values (e.g. axial mass) or 
to the approximation assumed in the Marteau 
model where the terms of order jj (p is tar- 
get nucleon momentum) are omitted from the 
hadronic tensor. In the detailed comparison of 
several MC codes predictions which was pre- 
sented by Zeller during her Nulnt02 talk 0] even 
bigger differences were shown. 

The results for the scattering on the oxygen 



nucleus are shown in Fig. 2. The NUX+FLUKA 
events are classified based on the particles leav- 
ing the primary vertex and as a result the FSI is 
not taken into account. To be more specific for 
the plot of the CC quasi-elastic scattering cross 
section only events with proton and muon pro- 
duced in the primary vertex are considered. On 
the other hand the plot for the Marteau model 
was obtained without the RPA corrections. In 
other words two free Fermi gas models are com- 
pared. 

The difference between the plots seen in Fig. 2 
has at least two reasons. Firstly it is inherited 
from the difference seen in Fig. 1. Secondly it can 
be caused by different distributions of the target 
nucleons momenta assumed in the two models. In 
the Marteau model a quadratic distribution with 
sharp cut at fcj? = 225 MeV is assumed while 
in the NUX+FLUKA event generator a smooth 
distribution is used |10j . 

To be able to compare theoretical predictions 
with experiment such nuclear effects as the rein- 
teraction or the A NN channel must be taken 
into account. 



3 



CD 
CO 

w 

CO 

o 

6 



12 



10 



Marteau model 
Nux+Fluka 

non-it decays of A 
Nux+Fluka - single pions absorption 

M » * f ♦ + + + + + + 



0.4 0.6 0.8 1 1.2 1.4 1.6 
Neutrino energy (GeV) 



1.8 2 



Figure 3. The cross section for the CC quasi- 
elastic scattering on 16 O (no pions in the final 
state). The contributions from pion absorption 
(NUX+FLUKA) and non-pion A resonance de- 
cays (Marteau model) are shown separately with 
the dotted lines. 



rises with the neutrino energy reaching about 30% 
at 2GeV and causes the cross-section produced by 
the NUX+FLUKA generator to be slightly higher 
than the one obtained in the Marteau model. 

3. SINGLE PION PRODUCTION 
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The plots in Fig. 3 were obtained by using the 
particles seen in the final state to select the CC 
quasi-elastic scattering events. This means that 
in the NUX+FLUKA generated event the pri- 
mary vertex might be classified as either quasi- 
elastic or single pion production but not only. 
However, we neglected rather small contributions 
from other primary vertex kinds because they are 
not included in the Marteau model. The plot ob- 
tained for the Marteau-Model includes RPA cor- 
rections. 

The comparison of plots found in figures 2 and 
3 leads to the conclusion that nuclear effects in- 
crease the cross section. 

This increase comes from an extra contribu- 
tion from pion absorption in the case of the 
NUX+FLUKA generator and from the non-pion 
decays of the A resonance in the case of the 
Marteau model. It is interesting that two con- 
tributions of so different origin are of the same 
order of magnitude. 

The difference between these two corrections 



Figure 4. Total cross section for single 7r + pro- 
duction on free protons. Experimental points are 
taken from [TT ] .[T2 ] . 



Predictions for single 7r production on free nu- 
cleons are compared only for the vp — > /i~p7r + 
channel. In this channel most of pions are pro- 
duced through A ++ excitation. In the experi- 
mentally measured hadronic invariant mass dis- 
tribution a resonance peak is clearly seen 
In the case of vn — » n~nir + and vn — > /i - p7r° 
channels the resonance peak is smeared out which 
indicates that the dynamics is more complicated. 

In Fig. 4 we compare the total cross sections 
predicted by both models with the experimental 
points. While both curves agree roughly with the 
results of Barish^T], the data points from J2| 
seem to favour the Marteau model. 

In Fig. 5 a similar comparison is made for the 
normalized hadronic invariant mass distributions. 
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Figure 5. Normalized invariant hadronic mass 
distribution in single 7r + production on free pro- 
tons. See explanation in the text. 



Figure 6. Total cross section for single 7r produc- 
tion on 16 (9 before Final State Interactions. 



The neutrino energy has been fixed at 1 GeV 
in both models because the experimental points 
were taken from |1 lj where most of the events 
came from neutrinos of energy around lGeV. In 
order to compare only distribution's shapes all 
the plots have been normalized to yield the same 
total cross section value (that of the Marteau 
model at neutrino energy lGeV). The shape ob- 
tained in the Marteau model agrees well with the 
data points which further justifies the hypothe- 
sis that most of the pions come from the A ++ 
decay. The disagreement between the shape ob- 
tained from NUX+FLUKA event generator and 
the data points is not surprising in this context as 
it is known not to contain any explicit resonance 
contribution. 

In the next two figures the total cross sec- 
tions of the pion production on oxygen are pre- 
sented. Contributions from all three pion pro- 
duction channels were added here because in the 
Marteau model they are related to each other by 
the isospin Clebsch-Gordan rules and they can- 
not be tested independently. In Fig. 6 we classi- 
fied the events based on the particles leaving the 
primary vertex. There is a major disagreement 



between the two curves: NUX+FLUKA predicts 
cross section to be much higher then the one ob- 
tained from the Marteau model. 

At neutrino energy 2GeV the difference is al- 
ready by a factor of 2 and it is seen that in 
the Marteau model the cross section saturates at 
much lower energy. The two main reasons for this 
discrepancy are the following: 

1) in the Marteau model pions are produced 
only via the A excitation with other (non- 
resonant, higher resonances) contributions ne- 
glected. Therefore the cross sections are underes- 
timated 

2) as shown in 4 the NUX+FLUKA predic- 
tions in vn — ► /i _ ?i7r + and vn — > /i~p7r° channels 
are significantly overestimated. 

In Fig. 7 the cross sections for single pion ap- 
pearance in the final state are shown. Prediction 
of the Marteau model is again much lower. Anal- 
ysis of 0] shows that already at neutrino energy 
2GeU NUX+FLUKA predict cross section higher 
then saturation plateau of approaches based on 
Rcin-Schgal model. By adding three contribu- 
tions one can expect that the plateau should be 
at approximate value 10 • 10 -38 cm 2 . 
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Figure 7. Total cross section for single 7r produc- 
tion on 16 O (Final State Interactions included). 



Figure 8. Fraction of pions produced from 
quasielastic primary vertex due to nuclear effects. 



Finally we considered the fraction of pions that 
are produced due to the Final State Interactions. 
In Fig. 8 we show the predictions of both ap- 
proaches. The value of this fraction obtained in 
the NUX+FLUKA scheme is four times larger 
then the one calculated in the Marteau model. 

This big difference can be discouraging but it is 
instructive that questions like this can be posed 
and answered within both numerical and analyt- 
ical frameworks. 

We hope that the discussion about numerical 
and analytical approaches to the description of 
nuclear effects in neutrino interactions can lead 
to improvements in existing MC codes. 
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